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Abstract 



A simple model of eddy currents in which current is computed solely from magnetic 
forces acting on electrons proves accessible to introductory students and gives a good 
qualitative account of eddy current forces. However, this model cannot be complete; 
it ignores the electric fields that drive current outside regions of significant magnetic 
field. In this paper we show how to extend the model to obtain a boundary value 
problem for current density. Solution of this problem in polar coordinates shows that 
the electric field significantly affects the quantitative results and presents an exercise 
suitable for upper division students. We apply elliptic cylindrical coordinates to 
generalize the result and offer an exercise useful for teaching graduate students how 
to use non-standard coordinate systems. 



1. Introduction 

Every student of Electricity and Magnetism learns that Lenz's Law predicts a force that 
opposes the motion of a conductor passing through a non-uniform magnetic field. Motion of 
the conductor's free charge through the field results in magnetic forces that drive current 
in the conductor. This current, in turn, interacts with the B field and results in a net 
magnetic force acting on the conductor. The current is called an eddy current. 

A classic classroom demonstration of eddy currents is a swinging metallic pendulum 
that passes through the field of a strong magnet. Eddy currents within the conductor 
damp the oscillation rapidly. When the conductor is replaced by another with holes, the 
eddy currents are impeded from circulating and the damping effect becomes very small. 
The currents cease in this case because a Hall electric field develops that balances the 
magnetic force acting on the free charge. 

We can estimate the eddy-current force acting on the conductor by using a few sim- 
plifying assumptions 1 ' 2 ' 3 . First, model the conductor as a very large plane sheet passing 
between circular magnet poles of radius a. Then, idealize the magnetic field as uniform 
in the cylindrical volume between the magnet poles and dropping abruptly to zero outside 
that volume. Figure 1 illustrates the model in a view perpendicular to the conducting 
sheet. 
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Figure 1: Magnetic field penetrating a circular portion of a moving conductor. 

In this view, the magnetic field falls on a circular region of the conductor, is uniform 
within the circle, and is zero outside the circle. The conductor moves in the x direction 

1 Susan M. Lea, John R. Burke Physics: The Nature of Things (Brooks/Cole, 1997), p. 974. 

2 R. K. Wangsncss, Electromagnetic Fields, (Wiley, New York, 1986) 2nd ed. See prob 17-14 p. 282. 

3 For a careful discussion of the much more sophisticated theory due to Maxwell and an extensive list of 
references see W. M. Saslow, "Maxwell's theory of eddy currents in thin conducting sheets, and applications 
to electromagnetic shielding and MAGLEV" Am. J. Phys. 60, 693-711 (1992). 
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with speed v and the magnetic field is given by B = — B z within the circle of radius a. 
By Ohm's Law, the current density is proportional to the force that drives it: 

— # — * 

J = (TV x B = avB z, (la) 

where a is the conductivity of the metal conductor. Then the force acting on a volume 
element dr of the conductor within the field is: 

dF = J x B dr. 

The net force acting on the conductor is 

F = -<tv B$Vx, (lb) 

where V is the volume of the conductor exposed to the field. 

This calculation correctly illustrates Lenz's Law and the dependence of the force on 
velocity and magnetic field strength. So, it gives a useful back-of-the-envelope estimate 
for the eddy current force. However, it is a somewhat naive estimate. Once the current 
leaves the vicinity of the B field, the model does not explain what causes the flow of free 
charge. It lacks an account of the Hall electric field (arising from charge distribution on 
the surface of discontinuity of B) which drives the current outside the B-field region, thus 

— * — * 

completing the current loops. This E field also opposes the current flow within the B-field 
region, indicating that equation (lb) overestimates the force. In this paper we develop a 
method to account for this effect, and so to improve the estimate. 



2. The Exact Circle Problem 

Calculating the charge densities that give rise to electric fields driving current in con- 
ductors is notoriously difficult 4 , but is usually not necessary. Here we can develop the 
calculation of current density as a two dimensional boundary-value problem using polar 
coordinates in the rest frame of the magnets. We retain the simple model of the magnetic 
field from the introduction and, for now, model the plate as infinite in the dimensions 
perpendicular to the field. We also assume that the plate's speed is sufficiently small that 
we can model the current distribution as a quasi-steady state in the magnet frame. The 
resulting problem is a challenging but accessible problem for upper division E&M students. 
In classic form, we observe that the current density is derivable from a potential that sat- 
isfies Laplace's equation except at the magnetic field boundary, develop the appropriate 
boundary conditions and solve via expansion in eigenfunctions. 

Then, the current density throughout is determined by Ohm's Law: 

J = a (E + v x B ) . 

4 J.D. Jackson "Surface charges on circuit wires and resistors play three roles", Am. J. Phys. 64, 
855-870 (1996). 
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Taking the curl of both sides and using a vector identity for V x (v x B) yields: 

Vx J = — ct(v ■ V)B. (2.1) 

Here we used the fact that the conductor's velocity and B are constant vectors in the magnet 
frame. Inserting the assumed form for B and using polar coordinates for V, equation (|2.1|) 
can be rewritten as: 

V x J = — av B cos((f))S(p — a)z, (2.2) 

where p is the radial coordinate with origin at the center of the magnetic field region. Next, 
since no charge buildup is expected with time, the equation of continuity demands that J 
be divergence free: 

V-J = 0. (2.3) 

Since the current density is curl free except at p = a, it is the gradient of some scalar 
potential on each side of a: J > = — V$> and J< = — V$<, where the subscripts > and < 
refer respectively to regions outside and inside the boundary. Thus V $ = everywhere 
except at the magnetic field boundary, and we may proceed with standard methods for 
solving Laplace's equation. 

The boundary conditions on the components of J perpendicular and parallel to the 
boundary follow from the divergence and curl of current density. This is a standard calcu- 
lation, with the results 

J>, x (a + ) - J<,x(a-) = 0; 

^<,ll( a+ ) - J<,\\{ a ~) = -OrV B COs((fi). 

In terms of potential, the boundary conditions are: 

d p <S> > (p = a + )-d p <5> < (p = a-) = 0; (5a) 
d^ y (p = a + ) — d^ < (p = a') = aav B cos (0). (5b) 

The problem is now completely specified, and we proceed by expanding the potential in 
eigenf unctions of the Laplace operator in polar coordinates. 

$(x) = { $< = ^"°=° pH ^ An sin ( n< ^ + Bn cos (W0] P < a ( 25 ) 
1 $> = Y^n=o P~ n l A n sin(n0) + B n cos(ra0)] p > a. 



From the boundary condition for Jy (eqn I5b|) and the orthogonality of the trigonometric 
functions, we see that the n — 1 terms are the only non-zero terms in the sums. Further- 
more, only the A coefficients are non-zero. The boundary conditions on the potential now 
give two equations for the coefficients. We find: 



$< = — \ctvoBqp sin(0) p < a 
$> = +lav B ^ sin(0) p > a. 



We can take gradients to obtain the current density: 



Jfx) 



-V$< = ±av B y 



\av B ( - ) sin(0)p - cos(4>)4> 



p < a 
p > a. 



(2.6) 



We still find a uniform current density within the magnetic field region. The current outside 
the field region follows a classic dipole pattern. The corresponding electric field that drives 
current in the region p > a and opposes it in the region p < a is found from 



E= v x B 

\v Q B (sp) sin(0) p - cos(0) 



p < a 
p > a 



(2.7) 
(2.8) 



The charge density a c that gives rise to this field is localized at p = a and is found from 
the standard boundary condition: 



a c = e p ■ 



(2.9) 

= e Q v B sm{<P) (2.10) 

For a field Bq = IT and a plate speed of vq = lm/s, the charge density is of order 
a c ~ KT 11 C/m 2 . 

Comparison equation ()2.6|) with equation (fTaj) shows that the current, and hence the 
net force acting on the conductor is half that predicted by the naive model. Such a simple 
result, in contrast with the complex correction one might have expected, rasies the issue 
whether a correction factor of 1/2 is generally correct or specific to the circuclar field 
geometry. We investigate that question in the following sections. 



3. Elliptical Magnetic Field Region 

The result for a circular magnetic field geometry demonstrates that electric field has 
a significant effect on eddy current flow. We were intrigued whether the factor of 1/2 
reduction is a general result or special to the case of circular geometry. To investigate 
this question, we solved the problem of an elliptically shaped magnetic field region with 
eccentricity e, as shown in Figure 2. The method follows the same outline as the circle 
problem except that we expand the potential in elliptic cylindrical coordinates, defined 5 ' 6 
in terms of Cartesian by: 

x = h cosh(?7) cos(ijj) y = h sinh(^) sin(ifj) z = z. (3-1) 

5 Morse and Feshbach, Methods of Theoretical Physics (McGraw-Hill, New York, 1953), Vol. 1, p. 514. 
6 Moon and D.E. Spencer, Field Theory Handbook, (Springer Verlag, Berlin, 1961), pp. 17-19. 
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The unit vectors are given by 

fj sinter/) cos(/0) — xj) cosier/) sm(ip) „ r) cosh(?7) sin(^) + ip sinh^) cos(^) 
J cosh 2 (77) — cos 2 (V0 J cosh 2 (77) — cos 2 (V0 

(3.2) 

The constant h is the product of the semi-major axis a and eccentricity e of the elliptical 
magnetic field region. The boundary of the magnetic field is defined by the level curve 

V — Vo = cosh -1 ( - J . (3.3) 



77=770 




Figure 2: The magnetic field penetrates an elliptical portion of a moving conductor. The 
level curves of coordinate rj are ellipses with foci p and q, a distance 2h apart. The level 
curves of tp are semi-hyperbolae; the value of -0 is the polar angle of the asymptote. The 
boundary of the magnetic field is given by rj — r]o = cosh _1 (i). The semi-major axis of the 
boundary is a. 



In these coordinates, the vector expressions for the divergence and curl of J are un- 
changed: 

V-J = V x J = -cr(v • V) B(rj) 

where B = —B z within the field region and zero outside. Using equation (J3.2)) we may 
express (v • V)B in this coordinate system. We find: 



V x J 



avo sinh(7?) cos(^) 9B(r/) 
h[cosh 2 (r]) — cos 2 (■?/>)] di] 
(JVqBq sinh(?7) cos(xp) 



5(77- 770) z. 



/i[cosh 2 (?7) — cos 2 ( , 0)] 
As before, conditions on the curl and divergence of J lead to boundary conditions: 

J>Avo) - J<Avo) = 0, 

av B sinh(r/ ) cos(tp) 



J>.\\(Vo) ~ J <,\\(Vo) 



cosh 2 (r/ ) — cos 2 (^) 
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Once again we can make the argument that J is the gradient of some scalar potential $ 
in the regions separated by 77 = i] . That is, J< = — V$< for i] < i] and J> = — V$ > for 
77 > i] Q . Components of current which lie perpendicular and parallel to the curve bounding 
the B region are then given by: 

1 <9$ 



J ± = -V$ • fj ■■ 
J = - V$ • t/> 



h\J cosh 2 (77) — cos 2 (■?/>) ^ 
1 .9$ 

h\J cosh 2 (77) — cos 2 (ip) ^ 



The boundary conditions on J then give us the following boundary conditions on the 
potential: 

d^yii] = - d^^i] = 770 ) = o-hv B sinh(i] ) cos(^), 
d v <S> > (r ] = r ] +)-d v <S> < (r ] = r ]o ) = 0. 

The expansion of the potential in terms of eigenfunctions of the Laplace operator in this 
coordinate system is 7 : 



$(x) 



= J2n=o i E n sinh(nr/) + F n cosh(n^)] [A n sin(n^) + B n cos(n^)] 77 < i] 



. = E"=o [G n e nri + H n e-^} [C n sin(n^) + D n cos(t^)] V > Vo . 

(3.4) 

Once again, orthogonality of the trigonometric functions ensures that only n = 1 terms 
will be non-zero. In the limit 77 — > 00, we require $> to remain finite, so we have G± = 0. 
The limit 77 — > describes the portion of the x-axis between the two foci. Here the curl and 
divergence of J are both zero, so that J is continuous across the x-axis. Since 77 increases 
away from the x-axis on both sides, fj changes direction discontinuously across the x-axis. 
Thus, continuity of J ± implies that J v oc d ri & < changes sign across the x-axis. This fact 
requires that the cos(?/>) term in $< be zero since only sm(ifj) is discontinuous at y — 0. 
Now, if) is also discontinuous across the x-axis between the foci, so continuity of Jy requires 
that d. i p§ < oc cos(ip) either be discontinuous (which it isn't) or be zero. Thus d^Q < = as 
77 — > 0, which requires that F\ = 0. Then, the boundary conditions at 770 require 



$(x) 




- costC+I h hfa) Sillll (^) Sill W V < VO 

ahv B sinh(77 ) cosh^^e - '' sin(^) 77 > r] . 



We can take gradients to calculate the exact current density: 

cosh(»7o)+sinh(»;o) 

J(x) = W _. n • , , \ __■,_/._ \ ( f)siii(^)-V>cos(i/>) \ _ ._ (3-5) 



J> = av B sinh(7] ) cosh(7] )e 57 T'"^ "Ti ^ > ^o- 

ycosn (??)— cos^(i/>) 



7 P. Moon and D.E. Spencer, op. cit., 19. 
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Substituting for r] in terms of the B-field region's eccentricity feqn !3.3|) we find: 



J J< = WoBo^Ls y n < cosh-^l/e) 

m = J> = av B ^e~v ( ) r, > cosh^l/e). ^ 

\ e \ y cosh {rj) — cos 2 (ip) J 

Again we find uniform current density in the magnetic field region. The factor of 1/2 
reduction in current density found for the circular field turns out not to be general. It is 
replaced by the factor 



i + VT 



e 



2 



which, of course, has the limit 1/2 as e — > 0. The graph of this function is shown as Figure 3. 
Since the force acting on the sheet is F net = — 5 |J < |yx where V is the volume exposed 
to the magnetic field, the force also has the expected limit. It is much more intricate and 
much less crucial to establish that the expression for J> reduces to the circular results in 
the limit e — ► 0. The calculation is not given here, but a copy of it is available from the 
authors upon request. 



4. The Effect Of Finite Conductor Size 

Once these two calculations are set up for infinite plates, it is easy to estimate the 
correction for finite plate size. One changes the boundary condition from |J>| — > at 
infinity to vanishing of the radial component of J> at a finite radial coordinate. 

For the circular case we take the plate to have a finite radius R. A solution is only 
feasible for the time when the plate is centered on the magnetic field region, so the result 
offers only an order of magnitude estimate of the effects. 

In eqn (|2.5|) (the expansion of the potential) an extra term in that increases with p 
is necessary to match the new boundary condition. A straightforward calculation reveals: 

lav B ( 1 - -^r ) y p < a 



R 2 

j(x) = r ) x2 y , (4.i) 

(R 2 - p 2 ) sin(0)p - (R 2 + p 2 ) cos(0)0 R> p> a. 



One may quickly verify that these expressions have the correct limits (eqn I2.6|) in the 
infinite conductor case (R — > oo). The effect on the dipole current term is substantial near 
the boundary. Current in the magnetic field region is further reduced by the edge effects, 
but by an insubstantial amount, unless the distance from the center to the nearest edge is 
comparable to the radius of the field region. 

A similar calculation is possible for elliptic cylindrical coordinates with a border at 
rj = H (semi-major axis of the boundary is A = aecosh(H). As in circular geometry, we 
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augment the old conditions with the new condition that the elliptic-radial component of 
J> vanishes at the boundary: d v ^{H,^) = 0. We find: 

J<(x) = av B o y-^ (\/ 1 "(f) 2 " v/r ^) 71 < cosh ' 1 ^- 

Observing that y/l — e 2 / (1 + y/l — e 2 ) = \/l — e 2 (l — y/l — e 2 )/e 2 , we see that the effect on 

J< is to replace 1 with \J 1 — (^) 2 in the last factor of the expression. Again the correction 
is of order of the square of the ratio of magnet size to the plate dimension. 



5. Conclusion 



We have demonstrated that the electric field has a significant effect on the eddy-current 
force computed from a simple model. The model gives the magnitude of the force as 

F = av B 2 Vf, (5.1) 

where V is the volume of conductor exposed to the magnetic field B , a is the conductivity 
of the metal and v o is its speed relative to the source of the magnetic field. The factor / is 
the correction due to the electric field; 1/2 for an infinite metal plate and circular magnet 
poles and y/l — e 2 / (1 + y/\ — e 2 ) for an infinite metal plate with elliptical magnetic poles. 
The first result follows from a boundary value problem accessible to an upper division 

_Vi^_£l_ 
l + Vi - 
.5; — 

0,4; 

0,3; \ 
0,2; \ 

, 1 : \ 



0,2 0,4 0,6 0,3 1 

Figure 3: The reduction factor in the current density for an arbitrary ellipse as a function 
of the ellipse's eccentricity. 



student, while the second result requires boundary value techniques that would be good 
training for a graduate student. Corrections for finite plate size alter the result by terms 
of order (magnet size / plate size) 2 . 

A possible objection to this method is the need for assuming an abrupt edge to the 
magnetic field region. Burke and Lea have developed a method for treating a more realistic 
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model of the field 8 . In the limit of zero separation of the magnet poles they find / = 1/2. 
For a pole separation of one-tenth of the pole radius, they find / = 0.39266. 

In all cases, < / < 0.5, though 0.5 seems to be a robust approximation. 
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